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Abstract 

We present a measurement of time-dependent CP asymmetries and an updated determination of 
the CP-odd fraction in the decay — > D*~^D*~. The measurements are derived from a data 
sample of 88 x W^BB pairs collected by the BABAR detector at the PEP-II energy asymmetric B 
Factory at SLAC. All results are preliminary. 
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1 Introduction 



The symmetry for combined charge conjugation (C) and parity (P) transformations is violated 
in B decays. Measurements of CP asymmetries by the BABAR |]l| and BELLE collaborations 
established this effect and are compatible with the Standard Model expectation based on the current 
knowledge of the Cabibbo-Kobayashi-Maskawa |^ quark-mixing matrix elements. 

As a result of the interference between direct B decay, expected to be dominated by the tree 
decay diagram, and decay after flavor change, a CP-violating asymmetry is expected in the time 
evolution of the decay^ B^ — > D*^ D*~ , within the framework of the Standard Model Up to 
corrections due to theoretically uncertain penguin diagram contributions |^], this CP asymmetry 
is related to sin2/3 {j3 = arg [ — ]). Penguin-induced corrections are predicted to be 

small in models based on the factorization approximation and heavy quark symmetry; an effect of 
about 2% is predicted by Ref. A comparison of measurements of sin2/3 from charmonium- 
containing b ccs modes such as B^ J/ijjKg, with that obtained in B^ D*~^D*~ is an 
important test of these models and the consistency of the Standard Model. 

The B^ — > D*~^D*~ mode is a pseudoscalar decay to a vector-vector final state, with contribu- 
tions from three partial waves with different CP parities: even for the S- and D-waves, odd for the 
P-wave. In the model described in Ref. |p the P-wave contribution is predicted to be about 11%. 
The angular distribution of the decay products can be used to measure the CP parameters of the 
CP-even and CP-odd components 

Following our initial results on this channel |1C], we present here an updated determination 
of the CP-odd fraction R±_ in the decay 5° D*+ D*- in BABM, based on a projected one- 
dimensional time-integrated angular analysis. We also present a preliminary measurement of the 
time-dependent CP asymmetry, obtained from a combined analysis of the time dependence of 
flavor-tagged decays and the one-dimensional angular distribution of decay products. 



2 The BaBAR Detector and Dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring. The data sample used for the time-dependent CP-asymmetry measurement corresponds to 
88.0 X 10^ e^e~ — > T{AS) BB events and the sample used for the measurement corresponds 
to 84.4 X 10^ -B-B pairs. The collider is operated with asymmetric beam energies, producing a boost 
(/97 = 0.55) of the T(4S') along the collision axis. 

BABAR is a solenoidal detector optimized for the asymmetric beam conflguration at PEP-II, 



and is described in detail elsewhere |11]. Charged particle (track) momenta are measured in a 
tracking system consisting of a 5-layer, double-sided, silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) filled with a gas mixture of helium and isobutane, both operating within 
a 1.5 T superconducting solenoidal magnet. Photon candidates are selected as local maxima of 
deposited energy in an electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals 
arranged in barrel and forward endcap subdetectors. In this analysis, tracks are identified as pions 
or kaons by the Cherenkov angle 9c, measured using a detector of internally reflected Cherenkov 
light (DIRC), and by the energy deposition, dP/dx, in the tracking system. The flux return of 
the magnet is instrumented with resistive plate chambers interspersed with iron (IFR) for the 
identification of muons and long-lived neutral hadrons. 

^Charge-conjugate modes are implied throughout this paper. 
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3 Event Selection 



mesons are exclusively reconstructed by combining two charged D* candidates reconstructed 
in a number of D* and D decay modes. Events are pre-selected by requiring that the normalized 
second Fox- Wolfram moment |12| of the event be less than 0.6. We also require that the cosine of 
the angle between the thrust axis of the reconstructed B and the thrust axis of the rest of the event 
be less than 0.9. These criteria are used to reject events coming from the e~^e~ — > cc continuum 
process. Backgrounds from u, d, s continuum processes are negligible in this analysis due to the 
presence of two charm particles in the final state. 

The and modes reconstructed are K^tt^, K^tt^tt^, ^ K- 

D° ^ i^°7r+7r-, D+ K-TT+TT+, D+ K°7r+ and D+ K-R+tt^. 

Charged kaon candidates are required to be inconsistent with the pion hypothesis, as inferred 
from the Cherenkov angle measured by the DIRC and the specific ionization measured by the 
SVT and DCH. No particle identification requirements are made for the kaon from the decay 

vr'^vr" candidates are required to have an invariant mass within 25 MeV/c^ of the nominal 
mass [^]. The angle between the flight direction and the momentum vector of the candidate 
is required to be less than 200 mrad, and the transverse flight distance from the primary event vertex 
must be greater than 2 mm. A mass-constrained fit is then applied to each surviving candidate, 
in order to improve the mass resolution of the K^tt'^ and K^tt~^-k~ combinations. 

Neutral pion candidates are formed from two photons in the EMC, each with energy above 
30 MeV; the invariant mass of the pair must be within 20 MeV/c^ of the nominal vr*^ mass, and 
their summed energy must be greater than 200 MeV. A mass-constrained fit is then applied to 
these t:^ candidates. The vr*^ from D*~^ — > D~^t:^ decay ("soft" vr*^), however, is required to have 
an invariant mass within 35 MeV/c^ of the nominal vr" mass and momentum in the T(4S') frame in 
the interval 70 < |p* | < 450 MeV/ c, with no requirement on the photon energy sum. 

and D+ meson candidates are required to have an invariant mass within 20 MeV/c^ of the 
nominal or mass. The same interval is used for all modes except K~'n~^TT^, which has 
a looser requirement of 35 MeV/c'^ due to the momentum resolution of the vr'^. 

D*^ mesons are reconstructed in their decays D*~^ D^tt~^ and D*+ D+vr*^. We include 
D*+D*- combinations decaying to {D'^tt+ , D^tt-) or (i:>°7r+, D^tt^), but not {D+ir^ , D~tt'^) due 
to the smaller branching fraction and larger expected backgrounds. The and D'^ candidates 
are subjected to a mass-constrained fit and then combined with soft pion candidates. A vertex fit 
is performed that includes the position of the beam spot to improve the angular resolution of the 
soft pion. 

If an event contains both a D*^ and a D*^ candidate, each is subjected to a mass constraint 
fit, and then combined to form a B candidate. B^ candidates with well reconstructed D* and D 
mesons are selected by means of a mass likelihood variable, jC-Mass, that includes all measured D* 
and D masses. This variable is defined in Eq. [l|: 
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'-'Mass ) X 

fcore G{AmD*+; AmD*p^^,aAm,orJ 

+ (1 - /core) G(Am£,.+ ; Am^,. 
fcore G{Amjj,-; Ami). 



+ 



{I- fcore) G{Amij,-; AmD*p^^,crAmt^u) 



(1) 



where G{x] fi, a) is a normalized Gaussian defined with mean = fj, and RMS = a; the subscript PDG 
refers to the nominal value |13]. For (Jmo '^^ use errors calculated candidate-by-candidate. The 
parameter fcore is the ratio of areas for the core and tail Gaussians. This along with (TAmcore 
CAmtaii determined from fitting the Am distributions in simulated signal events. The value of 
— ln(£ A/ass) is used to select signal in preference to background, with a different requirement used 
for each D decay mode combination. In an event where more than one B candidate is reconstructed, 
the candidate with the lowest —In(CMass) value is chosen. 

The primary variables used to distinguish signal from background are the energy-substituted 
mass. 



rriES 



Beam 



and the difference of the B candidate energy from beam energy, 

AE = E%- E%^^^ 

where the star indicates variables evaluated in the T{AS) center-of-mass frame. The fits described 
in the following sections are performed on events required to have |A£'| < 25MeV and tties > 
5.2GeV/c^ . The requirements on — \n.{CMass) and AE' were chosen to optimize S'^/{S + B), 
where S is the expected number of signal events and B is the expected number of background 
events. The optimization process was done with samples of simulated signal events and with 
generic BB and cc Monte Carlo simulated events. 

Figure |^ shows the events in the m^s projection after all selection criteria have been applied. 
The fit to this distribution uses a Gaussian and an ARGUS function |14] parameterization for signal 
and background, respectively. For the data sample corresponding to 84.4 x 10^ BB pairs, the fitted 
signal yield is 126 it 13 D*^D*~ events. 



4 Measurement of the CP-odd Fraction in D*+D*- 

In this section we present a one-dimensional angular analysis to determine the fraction, R±, of the 
P-wave CP-odd component in the vector-vector final state of the B^ — > 0*^0*" decay. In the 
transversity basis the following three angles involving decay products are defined (see Fig. |2|): 

• the polar angle 9i between the momentum of the it~ in the D*~ rest frame, and the direction 
of flight of the D*~ in the B rest frame, 

• the polar angle Otr between the normal, z, to the D*^ decay plane and the vr^ line of flight 
in the D*~^ rest frame, and 

• the corresponding azimuthal angle 4>tT- 
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Figure 1: The uies projection of the data for — > D*~^D*~ . These events are required to have 
lAE'l < 25MeV. The solid (dashed) hue represents the result (background component) of the fit 
described in the text. 




Figure 2: Representation of the "Transversity Frame" for the decay B D*^D*^ . The momenta 
of the D*^ decay products are represented in the B^ rest frame, while the momenta of the D*~^ 
decay products are represented in the Z?*"*" rest frame. The x direction is defined by the direction 
of flight of the D*~^ in the B^ rest frame. The {x, y) plane is defined by the momenta of the D*~ 
decay products in the B^ rest frame. The "transversity" axis, z, is orthogonal to the (x, y) plane. 
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The time-dependent angular distribution of decay products in the transversity frame for the 
mode 5° D*+D*- is given by @0: 

1 d^r 9 1 



'tr 



rdcos6'idcos6'trd0trdt 327r |AoP + l^yP + |^_l|2 

{4|^o|^ cos^ 6*1 sin^ Otr cos^ . 
+2|^ll P sin^ 01 sin^ 6*^ sin^ ' 
+2|^^psin^ 6*1 cos^ 6*^ 
+\/2Re(A|^o) sin 26*1 sin^ 6*^ sin 20tr 
— \/2Im(^Q^_L) sin 20i sin 20tr cos 0tr 
-2Im(^|^_L) sin^ Oi sin 26^^ sin 0tr} • (2) 

where ^0)^|ji^± Eire the time-dependent decay amphtudes in the transversity basis. For the 
decay A± = —A±, under the assumption there is no direct CP violation. 

In principle, 11 unknown parameters can be extracted from a full angular, time, and i?-flavor 
dependent analysis. These parameters are: Mi = \Ai{t = 0)|, 2 relative phases, and the complex 
CP parameters Aj = {q/p){Ai/Ai), with i =||,0, _L, resulting from the interference of mixing {q/p) 
and decay amplitudes {Ai). Considering the size of the current data sample, a simplified strategy 
is adopted to extract only the CP-odd component from the cos 9tr distribution. 

Ignoring the B flavor and integrating over time, 6i, and (ptr results in the one-dimensional 
differential decay rate: 

^ T^V = - ^i^' + ^os' ^tr (3) 
r dcosfc'tr 4 2 

with: 

^ Ml 

^ M2 + Mf + Ml ■ ^ ' 

The measurement of R± is based on an unbinned maximum likelihood fit of the cos 9tr distribu- 
tion, with a simultaneous fit to the mEs distribution. The probability density function for the rriEs 
distribution is given by the sum of normalized ARGUS |14| and Gaussian functions; the relative 
weight of each function is given by a signal fraction f^ig, which is allowed to float in the likelihood 
fit. The likelihood is defined as: 

^ = [fsig X ^mES ("^ES ,i ; CTruES ) ^ ^R± (^OS Oi] R±) + 

i=l,n (5) 

(1 - /sig)x-^mEs("^ES,i;K) X Hkg{cos9i;b2) 

where n is the number of selected events in the tties distribution, ^tues signal Gaussian for 

the mES distribution and ^mEs background ARGUS function with shape parameter k. J^r^ 

^Eq.5.44 in the reference appears with wrong signs in the last two terms and has been corrected in Eq. |^ of this 
paper. 
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refers to the probability density function for signal events, given by Eq. ^. The background shape 
is modeled by a polynomial in cos^tr^ 

^kgicOS Otr-, 62) = X (1 + 62 COS^ 0^,), (6) 

where N is the normalization factor. 

We categorize our events in three types: D*^D*^ {D^tt~^ , D^tt^), {D^tt^ , D^tt^), and 
{D~^Tr^ , D^7r~). We distinguish these three types of events because events with a neutral slow 
pion and events with a charged slow pion have different background levels and cos 9tr resolution 
(note: ^tr is the angle between the slow pion from the D*^ and the transversity axis). 

We allow different signal fractions for each event type. Thus, the parameters floating in the 
likelihood fit are: 3 signal fractions, the background parameter 62, 3 m-Es parameters (cr and mean 
of Gaussian fit, and k from ARGUS shape), and R_l. 

A fit to the dataset yields a value of i?_L = 0.096 ± O.OGO{stat), neglecting possible biases from 
angular resolution in 9tr and detector acceptance. Figure ^ shows the distribution of cos6'tr for 
candidates with mss in the range m-Es > 5.27 GeV/c^, with the result of the fit projected in the 
same region. The background component of the pdf is shown as the dotted line. 




cos(ej 

Figure 3: Likelihood fit result to the cos^tr distribution of the D*^D*^ events. The data points 
shown are from the region uies > 5.27 GeV/c^ and the solid line is the projection of the fit result 
in the same region. The dotted line represents the component of the pdf for the background. 

The experimental resolution of Btr biases the measurement of Rj_; the simulated distribution of 
residuals on 9tr is seen to have significant tails caused by mis-reconstructed events. The presence of 
these tails distorts the cos^tr distribution, and thus produces a bias on R±. A smaller bias on R± 
is produced by the transverse momentum dependence of the detector reconstruction efficiency. The 
presence of the two slow pions from the D* decays makes this analysis susceptible to this effect. 
The net bias from these two effects was estimated from a study on simulated data, and is found to 
be +0.028. This was determined using simulated events with a generated value of R± very similar 
to the value measured in data. The corresponding systematic error incurred for this correction is 
taken conservatively to be the full size of the correction, namely 0.028. The next largest systematic 
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uncertainty affecting tlie Rj_ measurement arises from the background parameterization (0.005). 
The total systematic uncertainty on R± is determined to be 0.03, giving as a prehminary result 
the final corrected value: 

R± = 0.07 ± 0.06{stat) ± Om{syst). (7) 

This is an update to a previous BABAR measurement pO| ], with a factor of three reduction in the 
statistical error {R± dependent). 

5 Measurement of the Time-Dependent CP Asymmetry in 

^*+^*- 

In this section we present a determination of the time-dependent CP asymmetry in — > D*~^D*~ . 
We perform a combined analysis of the cos 6tv distribution of all selected events and the time 
dependence of those with a flavor tag. The inclusion of the angular dependence enables us to fit 
for the CP asymmetries of the CP-even and CP-odd components separately. 

Although factorization models predict a small penguin contamination of the weak phase differ- 
ence Im(A/) = — sin2/3 Q, a sizeable penguin diagram contribution cannot a priori be excluded. 

Thus, the value of A/ = p4(7y different for the three transversity amplitudes because of 

possible different penguin-to-tree ratios. These contributions are explicitly included in the param- 
eterization of the flavor-tagged decay rates described here. 

The decay rate, /-|_(/_), for a neutral B meson tagged can be obtained from Eq. |^ 

as: 

-\At\/TgO 

/±(At) = — \0{l-^AV)^V[Ssm{AmdAt) + Ccos{AmdAt)]\, (8) 

where At = t^ec — hag is the difference between the proper decay time of the reconstructed B meson 
(i^rec) and the proper decay time of the tagging B meson (Stag); t^o is the B^ lifetime, and Am^ 
is the mass difference determined from the B^-B^ oscillation frequency. The dilution factor, V, 
describes the effect of incorrect tags, with AD accounting for possible differences in the mis-tag 
probability between B^ and B^. The O, C and S coefficients are defined as: 

O = ^[(l-i?x)sin2^tr + 2i?±COs2 0tr] 

C = ^[(1 -R±) ]~ 1^'' sin^ + 2r J~ I^^I' COS^ ^tr] 
4 i -1- |A+|^ i + |A_l|^ 

S = -l[il-R.)^^siu^e,.-2R,^f^cos^9,A. (9) 

These coefficients contain the explicit dependence on the transversity angle 9tr defined in the 
previous section, which provides a separation between the CP-odd (cos^ ^tr) and CP-even (sin^ ^tr) 
components. 

The CP-even parameters, |A+| and Im(A+), are related to Ay and Aq by: 
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Table 1: Efficiencies e^, average mis-tag fractions Wi, mis-tag fraction differences Awi = Wi{B^) — 
Wi{B^), and Q extracted for each tagging category i by using the i?flav sample. 



Category 


e{%) 


w (%) 


Aw (%) 


Q{%) 


Lepton 
Kaon I 
Kaon II 
Inclusive 


9.1 ±0.2 

16.7 ±0.2 

19.8 ±0.3 
20.0 ±0.3 


3.3 ±0.6 
10.0 ±0.7 
20.9 ±0.8 
31.5 ±0.9 


-1.5 ± 1.1 
-1.3 ± 1.1 
-4.4 ± 1.2 
-2.4 ± 1.3 


7.9 ±0.3 
10.7 ±0.4 
6.7 ±0.4 
2.7 ±0.3 


All 


65.6 ±0.5 






28.1 ±0.7 



1±|A+|2 MJ±M2 

Im(A||) 2 , Im(Ao) , ,2 

Im(A+) _ T+lx^^W + T^A^^ ^^^^ 



1±|A+|2 M2±M2 

It should be noted that this formulation does not take into account detector acceptance. Therefore, 
a combined fit to the At and cos ^tr dependence of the data will give an "effective" value of R± , which 
is not necessarily identical to the acceptance-corrected value from the time-integrated measurement. 

A measurement of CP asymmetries requires a determination of the experimental At resolution 
and the fraction of events in which the tag assignment is incorrect. A mis-tag fraction w reduces the 
observed CP asymmetry by a factor D = 1 — 2w. The mis-tag fractions and At resolution functions 
are determined from a sample, -Bflavj of neutral B decays to flavor eigenstates {D^*^^h^{h'^ = 
7r"'",p"'", and af) and J/t/j K*^{K*^ — > K~^tt~)) as for the sin2/3 measurement using charmonium 
decays, described in detail elsewhere 0. 

We use multivariate algorithms to identify signatures of B decays that determine the flavor of 
i?tag- Primary leptons from semileptonic B decays are selected from identified electrons, muons, 
and isolated energetic tracks. We use the charges of the best kaon candidates to define a kaon tag. 
Soft pions from D*'^ decays are selected on the basis of their momentum and direction with respect 
to the thrust axis of Stag- A neural network, which combines the outputs of these algorithms, takes 
into account correlations between different sources of flavor information and provides an estimate 
of the mis-tag probability for each event. 

Using the outputs of the algorithms and the estimated mis-tag probability, each event is assigned 
to one of four hierarchical, mutually exclusive tagging categories. The Lepton category contains 
events with an identified lepton, and a supporting kaon tag if present. Events with a kaon candidate 
and soft pion with opposite charge and similar flight direction are assigned to the Kaon I category. 
Events with only a kaon tag are assigned to the Kaon I or Kaon II category depending on the 
estimated mis-tag probability. The Kaon II category also contains the remaining events with a soft 
pion. All other events are assigned to the Inclusive category or excluded from further analysis 
based on the estimated mis-tag probability. The tagging efficiencies £i for the four tagging categories 
are measured from data and summarized in Table [l|. The figure of merit for tagging is the effective 
tagging efficiency Q = X]jei(l — '2'Wif'. This algorithm improves Q by about 7% (relative) over the 
algorithm used in Ref. Ill 
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The algorithm for vertex reconstruction and the determination of At are described in detail in 



Ref. [16| . The time interval At between the two B decays is calculated from the measured separation 
Az between the decay vertex of the reconstructed B meson and the vertex of the flavor-tagging B 
meson along the collision axis. We determine the z position of the -Brec vertex from the charged 
tracks that constitute the Bj-^c candidate. The decay vertex of the i?tag is determined by fitting the 
tracks not belonging to the i?rec candidate to a common vertex. An additional constraint on the 
tagging vertex comes from a pseudotrack computed from the i?rec vertex and three-momentum, the 
beam-spot and the T(4S') momentum. Events with a At error of less than 2.5 ps, and a measured 
|At| < 20 ps are accepted. 

We determine the parameters Im(A+) and jA+j with a simultaneous unbinned maximum likeli- 
hood fit to the At distributions of the i^rcc and -Bfiav tagged samples (Fig. ^. The At distribution 
of the -Bfiav sample evolves according to the known frequency for flavor oscillations in neutral B 
mesons. The observed amplitudes for the CP asymmetry in the i?rec sample and for flavor oscil- 
lation in the -Bfiav sample are reduced by the same factor (1 — 2w) due to flavor mis-tags. The 
At distributions for the i?rec and -Bfiav samples are both convolved with a common At resolution 
function. Events are assigned signal and background probabilities based on their ttt-es values. Back- 
grounds are incorporated with an empirical description of their At evolution, containing prompt 
(zero lifetime) and non-prompt components convolved with a separate resolution function |[l6[| . 

A total of 38 parameters are varied in the fit, including the values of Im(A+) and jA+j (2), the 
effective CP-odd fraction (1), the average mis-tag fraction w and the difference Aw between B^ 
and B^ mis-tags for each tagging category (8), parameters for the signal At resolution (9), and 
parameters for the background time dependence (7), At resolution (3), and mis-tag fractions (8). 
Because the CP-odd fraction is small, the parameters |Aj_| and Im(Aj_) are poorly determined. 
Therefore they are fixed in the fit to 1.0 and —0.741 |^] respectively. These are the values expected 
if direct CP violation and contributions from penguin diagrams are neglected. The changes in the 
fitted values of Im(A4.) and |A+| for different input values of Im(A_L) (varied between —1.0 and 
1.0) and |A_|_| (varied between 0.7 and 1.3) are taken into account as systematic uncertainties. The 
preliminary results obtained from the fit (Fig. ^ are as follows: 

Im(A+) = 0.31 ±0.43(stat) ±0.13(syst) (11) 
|A+| = 0.98±0.25(stat) ±0.09(syst). (12) 

The dominant source of systematic uncertainty comes from the variation on the value of A^ 
(0.09 and 0.02 respectively for Im(A+) and |A+|). Other relevant sources are the angular acceptance 
and resolution of the detector (0.06 and 0.08), the level, composition, and CP asymmetry of the 
background (0.07 and 0.02), the uncertainty on the SVT internal alignment and boost (0.03 and 
0.02), limited Monte Carlo simulation statistics (0.02 and 0.01), and possible differences between 
-Bflav and Bj-^c mis-tag fractions and resolution function parameters (0.01 and 0.01). The total 
systematic error is 0.13 for Im(A+) and 0.09 for |A+|. 

If the B D*^D*^ transition proceeds only through the b ccd tree amplitude, we expect 
that Im(A+) = — sin 2/? and |A+| = 1. To test this hypothesis, we fix Im(A+) = —0.741 Q and 
|A+| = 1 and repeat the fit. The observed change in the likelihood corresponds to 2.7 standard 
deviations (statistical error only). More data is needed to establish whether there are significant 
contributions from other processes, in particular, penguin diagrams. 
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Figure 4: From top to bottom: number A'^^o of candidate events in the signal region with a tag, 
number A''^o of candidates with a tag, and the raw asymmetry {N^u — N^q ) / {N^o + A''^o ) , as 
functions of At. The sohd curves represent the result of the combined fit to the full sample. The 
shaded regions represent the background contributions. 
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6 Summary 



We have reported preliminary measurements of time-dependent CP asymmetries and a measure- 
ment of the CP-odd fraction for the decay D*'^D*~ . The measurement of R± represents an 
improvement in the statistical uncertainty by a factor of almost three compared to previous mea- 
surements. The time-dependent asymmetry measurements still have large statistical uncertainties. 
These should be reduced steadily in the coming years as BABAR accumulates additional data, thus 
allowing useful tests of the Standard Model. 
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